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I. INTRODUCTION 



Quantum effects in semiconductor nanostructures are 
widely exploited in optoelectronic devices such as light 
emitting diodes or lasers. With the increasing demand 
of renewable energy supply, high efficiency photovoltaic 
devices were propose d, w hich also make use of semicon- 
ductor nanostructurea^H. However, these devices are not 
based on the standard IIl-V direct band gap materi- 
als conventionally used in optoelectronics, but on sili- 
con, which is most common in solar cell applications. 
It was shown that even though the suppression of di- 
rect optical transitions is somewhat relaxed in low di- 
mensional silicon structures, it is the phonon assisted 
indirect processes that dominate the optical respons^, 
and thus a suitable theoretical description of the pho- 
tovoltaic behaviour needs to include these processes re- 
sponsible for energy and momentum transfer, in addition 
to the general requirements of a quantum photovoltaic 
device, which were addressed in 4J. As an extension of 
this previous work, the present paper presents a micro- 
scopic description of phonon-mediated optical transitions 
in indirect semiconductor nanostructures, based on the 
non-equilibrium Green's function (NEGF) formalism and 
thus compatible with advanced quantum transport the- 
ories. 

The paper is organized as follows. After the deriva- 
tion of the general expressions for optoelectronic rates 
within the non-equilibrium Green's function theory, these 
are applied to the case of a direct gap bulk semicon- 
ductor and compared to the standard Fermi-Golden-Rule 
(FGR) result. Next, the description of phonon-mediated 
indirect transitions within the NEGF formalism is dis- 
cussed and the resulting expression for the generation 
rate again compared to the FGR, both in an analytical 
approximation and a numerical computation for a simple 
effective mass model of bulk silicon. Finally, the theory 
is extended to thin films in order to model photocurrent 
generation in ultrathin indirect gap semiconductor p-i-n 
junctions, which form important building blocks of future 
nano-photovoltaic devices. 



II. NEGF THEORY 

In order to be able to describe both optical transi- 
tions and inelastic quantum transport in semiconductor 
nanostructures within the same microscopic picture, a 
theoretical framework based on the NEGF formalism was 
developed in fi^. Before considering the specific case of 
photogeneration in indirect semiconductors, the aspects 
of the theory related to interband transitions will be out- 
lined below in a general formulation, which relies on the 
extensive work on similar systems investigated with focus 
on their light-emitting and lasing propertie^^^'i^. 

A. Interband generation and recombination rates 

The macroscopic equation of motion for a photovoltaic 
system is the continuity equation for the charge carrier 
density 

atp,(r) + V-j,(r) =g,(r)-7^,(r), c^e,h (1) 

where pc and jc are charge carrier and particle current 
density, respectively, Qc the generation rate and TZc the 
recombination rate of carriers species Pl. The micro- 
scopic version of this equation in terms of NEGF corre- 
sponds tcPESl (1 = {n , G C}) 

lim [ih A + G(l,2) + [iJo(ri) - H^{v^)] G(l,2)} 

= lim / d3[E(l,3)G(3,2) - G(1,3)S(3,2)] , (2) 

where Hq is the Hamiltonian of the non-interacting 
electronic system, G is the electronic non-equilibrium 
Green's function defined on the Keldysh contour C and 
S is the self-energy encoding the interaction of the elec- 
tronic system with phonons, photons and with itself, i.e. 
the scattering processes which may give rise to intra- or 
interband transitions. In steady state, (jlj becomes 

V-je(r) -gc(r)-7ec(r), c = e,h. (3) 

In the microscopic theory, the divergence of the electron 
(particle) current corresponds to the limit 1 — 2 of the 
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RHS in (ip, 
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If the energy integration is taken over the range of all 
bands connected by the interband scattering process, i.e. 
both valence and conduction bands, the resulting diver- 
gence should vanish, if the self-energies are chosen prop- 
erly, ensuring the conservation of the overall current, 
which is nothing else than the corresponding formulation 
of the detailed balance requirement. If the integration is 
restricted to one of the bands (e.g. B=valence band or 
conduction band), the above equation corresponds to the 
microscopic version of (|3]), and provides on the RHS the 
total local interband scattering rate, 



Rrad{v) =n{v) - g{v) 

2 f dE 
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2Trh 



(5) 

S^(r,r';£;)G<(r',r;£;) 
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(6) 



Depending on the nature of the interaction described by 
S, the scattering process may be highly non-local, in 
which case the above rate contains contributions from 
a large volume. The total interband current is found by 
integrating the divergence over the active volume, and 
is equivalent to the total global transition rate and, via 
the Gauss theorem, to the difference of the interband cur- 
rents at the boundaries of the interacting region. Making 
use of the cyclic property of the trace, it can be expressed 
in the form 



Rrad I <rr 



B(r) 
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j:>{r,r';E)G<{r',r;E) 



T.<{r,r';E)G>{r',r;E) 
(7) 



with units [i?rad] = s ^. Within the general (basis- 
independent) NEGF picture. The RHS of Eq. ^ can 
be interpreted as follows: S~^I]<^^^(iJ) represents the 
rate at which charge carriers with energy E may leave 
(occupy) a state at that energy whereas G>'^<^E) quan- 
tifies the energy resolved probability that the system can 
accept (donate) a particle of energy E, i.e. that there 
is an unoccupied (occupied) state at the right energy. 
The first term thus represents the total inscattering rate, 
while the second term provides the total outscattering 
rate. 

If we are interested in the interband scattering rate, 
we can neglect in Eq. ([t]) the contributions to the self- 
energy from intraband scattering, e.g. via interaction 



with phonons, low energy photons (free carrier absorp- 
tion) or ionized impurities, since they cancel upon en- 
ergy integration over the band. However, if self-energies 
and Green's functions are determined self-consistently as 
they should in order to guarantee current conservation, 
the Green's functions are related to the scattering self- 
energies via the Dyson equation for the propagator, 

G«(^)(ri,ri,;i?) ^G^,^^\v,,vr; E) 

+ j d'r.J d\,G^^^\r,,r,;E) 

xS«(^)(r2,r3;£;)G^(-4)(r3,ri,;i?), 

(8) 

and the Keldysh equation for the correlation functions, 

G^{r,,ry;E)^ J d^^ j d\^G''{v^,V2;E) 

X E^(r2,r3;S)G^(r3,ri,;i;), (9) 

and will thus be modified due to the intraband scattering. 
This means that in the case of self-consistent solutions, 
it is in general not possible to completely separate the 
effects of the different scattering processes, nor to isolate 
coherent from incoherent transport. 

In the remainder of the paper, the general theory of in- 
terband transitions outlined above will be applied to op- 
tical interband transitions in direct and indirect semicon- 
ductors. For computational purposes and to ease com- 
parison with existing descriptions, the theory will be re- 
formulated using a simple effective mass band basis for 
completely homogeneous bulk systems and for inhomo- 
geneous thin film devices. 



B. Bulk semiconductor 

For a homogeneous bulk system, the field operators for 
carriers in band 6 can be written in the Bloch state basis. 



^b{r,t) = ^V'&k(r)cf,k(^)■ 



(10) 



The expression for the total radiative rate of carriers in 
band b is simplified by using the Fourier space represen- 
tation of Green's functions and self-energies, which reads 
(0 = G,E) 

0,,b. (r, r'-E)^Y. ^bi.ir)Of,^b' (k; £;) Vt* k(r'), (H) 



where Ot t,' (k; E) is the steady-state Fourier-transform of 
the real-time Keldysh components of the contour-ordered 
Green's function 

Ob,b'{k; t~t')^^ {Tc {^fak(i)4'k(i')}) ■ (12) 
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Inserting these expressions in ([7| , the band-resolved rates 
are obtained as 
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I]> (k;i?)G<,(k;i?) 



(13) 

In the following, the off-diagonal terms will be neglected 
(Ob = Ob,b'), which means that only incoherent inter- 
band and subband polarization is considered. 



1. Direct interband transitions 

Inserting the electron-photon self-energy for a two- 
band model of a direct semiconductor (App. |A]), and ne- 
glecting intraband processes (free-carrier absorption and 
emission), the expression for the interband absorption 
rate becomes (in the following, the energy integration is 
restricted to the conduction band) 

^-^=^/SEEl^^eX(k,A,q)|X. 



k A,q 



(14) 



with Al^y the optical matrix element (App. vKn and q 
the occupation of the photon modes. The Tatter is ob- 
tained from the modal photon flux via = (f^ qW^j 
where c is the speed of the light in the active medium. 
The modal photon flux in turn is given by the modal in- 
tensity of the EM fleld as (/)]| q = ^/ (/kj> q). Similarly, 
stimulated interband emission reads 

2 
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k A,q 

xG>{k;E-huj^)G<{k;E), (15) 
while spontaneous interband emission is expressed as 
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k A,q 

xG>{k;E-nuq)G<{k;E). (16) 

Rem,st of photons 



The net absorption Rabs.net = Rabs 
in mode (A, q) is thus given by 



i?a...„et(A,q) ^- j 2^El^^™(k,A,q)p7V,^_, 



G<{k-E-hLo^)G>{k-E) 
-G>{k;E~hLo^}G<{k;E) 



(17) 



Conventionally, the rates for absorption and emission 
are calculated based on Fermi's Golden Rule (FGR), cor- 
responding to the Born approximation within flrst order 



perturbation theory. It is thus instructive to compare 
the above expressions with the FGR-rate for net direct 
interband absorptiorP^^ 



-^abs.netCA, q) — 



27r 



k 



|Ar(k,A,q)|27V; 



A,q 



X S{ecik) - £,(k) - fic^q) [/„(k) - /,(k)] . 

(18) 

Here, the absorbing (bulk) material is described by the 
dispersion relations £b(k), 6 = c, w, and assumed to be 
in a quasi-equilibrium state with occupation described 
by the Fermi function fb{E) = (exp[/3(i? — fib)] + 1)~^, 
/3 = {ksT)"^, where /^b, b ~ c,v are global quasi-Fermi 
levels. 

In order to reproduce the FGR result from the more 
general expression in terms of Green's functions, the un- 
perturbed (i.e. non-interacting) equilibrium form of the 
latter needs to be used, which corresponds to the expres- 
sions for free fermions in equilibrium, 

i<(0) 



G<*'" (k; E) ^2mME)5{E - e^ik)) , 



G>("^ (k; E) =2m [ME) -l]S{E- eb(k)) , 



(19) 
(20) 

Introducing these expressions in ( |17[ ) and carrying-out 
the energy integration reproduces the FGR expression 

To estimate the deviation of the rate from the FGR 
result for self-energies beyond the flrst Born approxima- 
tion, the FGR rate is flrst used to derive the standard 
expression for the bulk absorption coefficient of the two 
band effective mass model, which amounts to 



a(/iw^) 



Rabs,net{f^t)/V 



(21) 



where Slhuj-y) is the monochromatic energy flux density 
of the EM fleld (i.e. the absolute value of the Poynting 
vector) given by 



where 



An ^ 

A 



27r2(;iS)3 



(22) 



(23) 



is the photonic density of states of an optically isotropic 
medium with refractive index rib — \/£b and correspond- 
ing speed of light c = co/n. With the standard approxi- 
mation of isotropic and momentum-independent optical 
matrix elements, i.e. |M^^(k, A, q)p « M.1^(hxi^), the 
absorption rate is rewritten as follows. 



Rabs{ti^l) ^M2y{huj^)Jcv{tuj^)py{hujj) 

dn 



47r ^ 

A 



(24) 
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which provides the bulk absorption coefficient 



(26) 



with = MJ^/c. The difference between FGR and 

NEGF approaches concerns the term Jcvifi^-y), which in 
the FGR case takes the specific form (using the contin- 
uum approximation ^j^. — ^^^.^3 J cPk ) 



n-FGR 
^ cv 



2tt 2 
'7r(27r)3 
X [/„(k)-/,(k)] 



y" d^fc (5(ec(k) - £„(k) - hhj^ 



(27) 



If the occupation depends only marginally on crystalline 
momentum, the above expressions related to the joint 
density of states Jcv with suitable occupation, 



jFGR 



27r 



(fi^^) =^Jcv{f^-,){f, - fc), (28) 

(29) 

In the NEGF case, joint density of states and occupa- 
tion cannot be separated, but are both contained in the 
Keldysh Green's functions 



2. Phonon assisted interband transitions 

Even in the case of a direct semiconductor discussed 
above, the presence of phonons can have a considerable 
effect on the optical transition rates as new excitation 
paths become available, i.e. the phonons strongly in- 
crease the number of initial-final state pairs separated by 
a given transition energy. However, comparing to the di- 
rect transition, this enhancement of the joint density of 
states is overcompensated by the fact that the transitions 
become much more unlikely due to the need for coupling 
to a suitable phonon and the detour via the virtual state. 
As a consequence, phonon assisted transitions may be 
neglected in direct bulk semiconductors where crystalline 
momentum is conserved. Obviously, the situation is com- 
pletely different in indirect bulk semiconductors, where 
no direct transitions are possible. There, at lowest non- 
vanishing order, four different excitation processes exist 
(Fig. [1]): 

5*1+ : 1. phonon absorption, 2. 

51- : 1. phonon emission, 2. 
5*2+ : 1. photon absorption, 2. 

52- ■ 1. photon absorption, 2. 



(27r)3 



dE 
2Trh 



(30) 
(31) 



which via Vcv,o = 'Pcv ~ ^™ related to the Keldysh 
components of the free-carrier interband polarization 
function 



(27r)3 



d-^k 



xGfik-q;E-E). 



(32) 



This expression is valid in any situation that can be de- 
scribed in terms of single-particle Green's functions. The 
deviations from the FGR result are marginal in the quasi- 
equilibrium case and in absence of further interactions 
beyond electron-light coupling and of sources of non- 
equilibrium, which both modify the Green's functions, 
causing them to differ from the expressions given in Eqs. 



(191 and (20). It is straightforward to show that insert- 



ing the latter expressions in (30) reproduces (27), and for 



the special case of spherical bands, the joint density of 
states is given by the well-known analytical expression 



) ri 

direct band gap. 
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-^E^-Eg 
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where rnl is the reduced effective mass and Eg is the 
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FIG. 1. Possible excitation pathways for optical interband 
transitions in indirect semiconductors on the example of bulk 
silicon. 



To exemplify the inclusion of the phonon scattering in 
the NEGF description of optical interband transitions, 
we will in the following focus on processes S2±- The 
FGR transition rate for these processes are obtained from 
second order perturbation theory as 

2^ ... ^ ^ [Mt^{K,\^.)?[Ml^{].,)Y 



jdFGR 
^abs,± 



(A,q) 



-N 



N- 



Zq X! X! 

ke,k„ A,Q 



1 1 

2 2 



ph 
A,Q 



|ec(k„) - et,(k„) - /iWqp 
/„(k„)[l-/e(k,)] 



X(5(ec(kc) - Si,(ky) ~ huj^ =f ^^^a,q)- (34) 
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"virtual" CB (r^) states 
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FIG. 2. Self-consistent computation of Green's functions and 
scattering self-energies in silicon enabling the description of 
phonon-assisted indirect optical transitions 



Here, M^'^q encodes the matrix element for coupling of 
electrons to the phonon mode (A, Q) with energy /if^A^q 
and occupation given by the Bose-Einstein distribution 



Q = (e'^''"'''^ - 1) , /3 = (ksT)-^ at lattice tem- 
perature T. In analogy to the direct case, with the pre- 
viously used approximations for the optical matrix ele- 
ments and the additional assumptions 



ISciK) ~ SviK) - fl^Cll ~ \EgQ - hUJc, 



(35) 



Ml'^^iK, ke) « Mf6{K - + Q) (36) 
we can write the (phonon-assisted) absorption coefficient 



\EgQ - fUxJ^l'^ 



ph 
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ph 
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2 2 



X^'ihiu,), (37) 



but where now the joint density of states for indirect 
transitions is used in (27 1, 



Jd^'h j d?k2 5(ec(ki)-£.(k2) 

(38) 



(27r)6 



which for spherical bands and a single phonon frequency 
may be simplified to 



-rind 
''cv,± 



(27r;i)3 ^ ^ 



Eg T hnf, 



(39) 



where Eg denotes the (indirect) band gap. 

In the following, the phonon-assisted absorption rate 
shall be derived within the NEGF formalism, start- 
ing from the expression for the absorption rate, but 
where now the Green's functions contain the contribu- 
tions due to the electron-phonon scattering. Since the 
Green's functions and interaction self-energies are eval- 
uated in a self-consistency iteration process, all possi- 
ble single phonon processes are included to all orders, 
i.e. the Green's functions contain contributions from a 
number of scattering processes that corresponds to the 
number of self-consistency iteration steps. It is thus 
via this self-consistent computation that phonon assisted 



optical transitions are enabled, even though the self- 
energies themselves are only on the level of a first order 
self-consistent Born approximation, i.e. do not include 
the combination of both electron-phonon and electron- 
photon scattering at the same time. In Fig. [2j this pro- 
cedure is shown for the technologically relevant exam- 
ple of indirect interband absorption of photons in silicon, 
where zone-boundary phonons provide the wave-vector 
difference in a F — X inter-valley scattering process. 

We start again from the general expression for the in- 
scattering rate, which in this case reads 



(40) 



Assuming a photon-first indirect process (5'2±), the cor- 
responding expression for F — X-scattering is inserted for 
the self-energy (see Eq. (A28) in App. A 2), leading to 



the equivalent of Eq. (14) {<j — LA, TO: phonon mode). 



R.n=l J ^^^[A/r(l].)]2[ArfG<(Q;i?-;.l]„ 



(TVf + 1)G< (Q; E + hn„)G>^ (k; E). 



(41) 



In the next step, the Keldysh equation for electron- 
photon interaction is used to replace the lesser F^-GF, 

G< (k;^;) = G^^(k;i;)E<J(k;£:)GJ^^(k;i;) (42) 

providing the modal absorption rate 

Rats{X,ci) = - J — ^^[Aff (f^.)]2[M^(k,A,q)]2 



k Q,(T 



K^G^^ (Q; E - hn^)G^^ (Q; E - hn„ - hoj^) 
X G^^(Q; E - hn„) + {NP^ + 1)G^^ (Q; E + hil,) 

X G^jQ;E + hn^ - hujq)G^^{Q] E + nn, 

xG>^{k;E). 



(43) 



Again, this is to be compared to the FGR result, which, 
for the same electron-photon and electron-phonon inter- 
action Hamiltonian terms, follows from second-order per- 
turbation theory to 



uFGR^, ^ 2^ ^ ^ [M^^f^.)]'[M^(k„,A,q)]Xq 
R-abs (A, q) =-r- 2^ 2^ 



kc,k„ 



|£r, (k„) - er„ (k„) - huj^ 



X [NfS{ex^ (k,) - £r„ (k„) - hco^ - hn,) 

+(A^f + l)S{ex^ (kj - er„ (k,) - hu^ + hn^ 
xfrMil-fxA^c)]. (44) 

Now, inserting the non-interacting equilibrium expres- 
sions for the lesser, greater and retarded GF in (43) pro- 
vides the expression 
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k Q,(T 



NP^ifr^ {E - nn^ - hujq)2TTS{E - hVl„ -hu:^- er„ (Q)) 



\E - hfla - - er,(Q) + iVQ+\ 



{NP'' + l)i/r„ {E + Ml„ - hujq)2TrS{E + hQ^ - hoj^ - £r„ (Q)) 



i[fxSE)-l]2nS{E~SxA'^)), (45) 



which, upon performing the energy integration and with 
Q = kt,, k = kc, reproduces again the FGR result given 



extension, the envefope basis function can be written as 



in (44) 



For the numerical implementation, a simple three-band 
effective mass model for the electronic structure of silicon 
is used. The X electrons are described by an multi-valley 
picture with identical values for transverse and longitudi- 
nal effective mass, and the holes as well as the virtual 
Fc states used in the indirect transitions are modelled 
by single effective mass bands. The holes are modelled 
by a single effective masses corresponding to heavy and 
light holes. The band structure and interaction parame- 
ters used in the numerical examples are listed in Tab. H] 
Fig. [3] displays the close agreement between NEGF and 
FGR for equilibrium bulk absorption close to the indi- 
rect band edge and the characteristic fingerprints of the 
involved TO and LA phonon modes. 



Thin film devices 



(r) 



Va 



(48) 



where ry = (x, y), A is the cross sectional area and Xi is 
the longitudinal envelope basis function. For the l atter, 
finite element shape functions are a popular choicd^J^i^. 
Here, we will use a simple finite difference basis equiva- 
lent to a separate single band tight-binding approach for 
each bancpHU], 



X^{z) = Mz - z,) - e{z - z,+i)] /VA. (49) 

In the above basis, the Green's functions and self-energies 
have the potentially non-local representation (O = G, S) 

Ofc,b, (r, v'-E)^Y.T. ^'fk,, {v)0,b.,w (k|i ; i?)^,Vk|| (r'), 

k|| i.J 

(50) 



In the case of ultra-thin-film and especially hetero- 
multi-layer devices, translational invariance does no 
longer apply in growth direction. For such a system, 
the appropriate representation of the field operators has 
the form 



k||,i 



with the basis functions 

■0mk||(r) = (/3ik|| (r)w„ko(r), 



(46) 



(47) 



where (piky is the envelope basis function for discrete spa- 
tial (layer) index i (longitudinal) and transverse momen- 
tum k|| , u„ko is the Bloch function of bulk band n, cen- 
tered on ko. In the case of a system with large transverse 



TABLE L Material parameters used in simulations 



■m*x 


mf„ 




Eg,Tv-X 


P'L/mo 


0.3 mo 0.3 mo 


0.54 mo 


3.1 eV 


1.17 eV 


4 eV 


a Mode 










(r-X)i LA 

(r-x)2 TO 


18.4 meV 2.45x10* eV/cm 
57.6 meV 0.8 xlO** cV/cm 







where the contour-ordered steady-state Green's functions 
are now defined as 



ib,jb' 



(k|| ; ^ - - ^ {Tc {^.k,, (t)c],,k„ it')}) ■ (51) 
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FIG. 3. Photon absorption rate for indirect optical transitions 
in intrinsic bulk silicon: FGR and NEGF rates are in close 
agreement, since the correction fo the conduction band DOS 
due to the phonon-mediated F — X scattering is weak. 
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FIG. 4. Convergence of tiie interband photogeneration rate 
(integrated over growtli dimension) and of the resulting to- 
tal interband current density in the NEGF self-consistency 
iteration process (Fig. [2]|. 



The conservation law corresponding to Eq. ([4| for the 
divergence of the electron charge current between model 
layers i — 1 and i becomes 



hAA ^ I 2-K 

k 



(52) 
(53) 



where Rrad{z) is again the local radiative rate. Making 
use of the cyclic property of the trace, the total radiative 
rate follows as 



Rrad 



^XfTr{i:[^e<(k;i.)G>(k;i.) 

-S>(k;i?)G<(k;i?)]} (54) 
^AY,Rrad{zi)=Y.{J^ -Jt,) (55) 



— jB tB _ J Jn, B — CB, 



Ji, B = VB. 



(56) 



In the case of the indirect gap material the total radiative 
rate is 

" k|| 

-S>_p^(kj|;i?)G< (k||;i?)]}, (57) 
and the inter-valley phonon scattering rate reads 



'ep(T-X)X 



2 f dE r 

k|| 

X G> (k|| ; E) - S> (k,, ; i?)G< (k,, -E)]]. 




.1 1.12 1.14 1.16 1.18 1.2 1.22 
photon energy Epj^^t [eV] 



FIG. 5. Photon absorption rate for indirect optical transi- 
tions, computed with the thin film model (dots) for a silicon 
p-i-n diode with low doping levels (10^® cm~^, in comparison 
with the Fermi Golden Rule rate for bulk silicon (full line). 



The convergence of scattering rate and interband current 
as a function of self-consistency iteration steps is shown 
in Fig. |4j The convergence is fast due to the low carrier 
density in the absorbing quasi-intrinsic region. 

Fig. [5] shows the indirect optical absorption for an 
quasi-intrinsic slab of the indirect gap semiconductor dis- 
cussed above. In this case of near flat-band conditions, 
the bulk absorption is retrieved close to the band edge. 
For the more interesting case of a ultra-thin bipolar junc- 
tion with high doping, the global absorption coefficient 
starts to deviate from the bulk value due to the strong 
internal field. This so-called Franz-Keldysh effect is dis- 
played in Fig. [6] Fig. [7] shows the energy resolved 
local generation rate in the virtual Fc-states and result- 
ing local photocurrcnt spectrum in the X-valley band 
at a photon energy -Ep^ot = 1-17 eV and very low in- 
ternal field F = 2 kV/cm. The generation is relatively 
uniform throughout the device, as well as the increase 
of electron and hole current components towards the re- 
spective contacts. Fig. |8] shows the corresponding sit- 
uation for the diode with strong internal field F = 169 
kV/cm. In this case, the photogeneration and current 
contributions are distributed over a large energy range, 
the latter due to the absence of efficient relaxation mech- 
anisms that would confine current flow closer to the band- 
edge. In the vicinity of the n-contact, the high electron 
concentration results in large electron-phonon intraband 
scattering current contributions (Fig. |9^), which due to 
the opposite signs of in- and outscattering components 
cancel upon energy integration over the band. The net 
current can thus be attributed exclusively to interband 
transitions, and the local sum of electron and hole con- 
tributions is perfectly conserved, as shown in Fig. |9b. 
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FIG. 6. (a) Band profile of a silicon p-i-n junction from the self-consistent solution of the NEGF-Poisson equations with the 
indicated doping levels from Nd ~ 10^® cm~'^ to Nd = 10^* cm~^. The corresponding fields in the intrinsic region vary from 
F = 2 kV/cm to _F = 169 kV/cm. The equilibrium Fermi level lies at 1.1 eV. (b) Near band gap optical absorption rate in the 
intrinsic region of a silicon p-i-n junction, for various values of the built-in field corresponding to different doping levels. 



Since the photogeneration rate is intimately related to 
the T — X intervalley electron-phonon scattering rate via 
the self-consistent computation of Green's functions and 
self-energies, and the photocurrent results from photo- 
generated excess charge that is transferred via the scat- 
tering process with phonons to the extended X-valley 
states, the rates of this phonon-mediated charge transfer 
are expected to show the signature of the photogenera- 
tion process. Indeed, if the rate of electron scattering out 
of (into) the virtual Fc-states into (out of) the X-states 
is considered, there is a large net scattering from Fc to X 
under illumination, while in- and outscattering rates are 
balanced in the dark. This is shown in Fig. [lO] for the 
case of the low-doping diode where the intrinsic region 
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FIG. 8. (Color online) Same as Fig. [7] but for a strong inter- 
nal field of F = 169 kV/cm. The large deviation close to the 
n-contact is due to the electron-phonon intraband scattering 
current that vanishes upon energy integration over the band 
(white=negative, black=large positive, see Fig. [9a). 
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FIG. 7. (Color online) Energy resolved local generation rate in 
the virtual F^-state and resulting local photocurrent spectrum 
in the X-valley band at very low internal field F = 2 kV/cm. 
The grey lines indicate the band edges (X- valley for electrons). 



FIG. 9. (Color online) (a) Current spectrum of X-valley elec- 
trons close to the n-contact. (b) Integrated net current for 
electrons (el) and holes (hi) in the high-field device. The to- 
tal current is perfectly conserved over the whole device. 
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FIG. 10. (Color online) Phonon mediated electron transfer rate between virtual (Fc) and extended (X) conduction band states 
(a) in the dark and (b) under illumination with Ephot ~ 1.17 eV, for the low-doping diode, where the intrinsic region extends 
over the whole device. While in-and out-scattering are balanced in the dark, there is a strong net charge transfer from Fc to 
X under illumination, and the inscattering rate reflects the spectral pattern of the photogeneration. 



'inscattering" 




38 40 42 
position z [nm] 



FIG. 11. (Color online) Phonon mediated electron transfer rate between virtual and extended conduction band states (a) in 
the dark and (b) under illumination with Ephot = 1.17 eV in the center of the short intrinsic region of the high-doping diode. 
Closer to the contacts, the scattering rate contributed by the photogeneration is masked by the large intraband relaxation rate. 



extends over the entire device. In the case of strong dop- 
ing, the electron-phonon scattering rate in the dark is 
much stronger close to the contacts, such that the effect 
becomes only visible in the central, intrinsic region of the 
junction. There, the phonon-related charge transfer rate 
again reflects the local photogeneration rate, as displayed 
in Fig. [n] 

III. CONCLUSIONS 

We have presented a novel theoretical approach to 
the description of phonon-mediated photogeneration in 
indirect semiconductors. While compatible with the 
standard Fcrmi-Golden-Rule approach in the quasi- 
equilibrium bulk limit, its range of validity extends to 
quantum transport in open systems involving arbitrary 



heterostructure states far from equilibrium and the ef- 
fects of non-locality in the scattering processes, which 
are important aspects of advanced photovoltaic and light 
emitting devices, where often spatial and spectral resolu- 
tion are required to gain access to a deeper understanding 
of the device characteristics. Thanks to this versatility, 
the theoretical framework lends itself to the modelling of 
indirect semiconductor based nanostructures with poten- 
tial applications in a multitude of optoelectronic devices, 
such as silicon-based quantum well solar cells. 
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Appendix A: Interaction self-energies 

The self-energies as defined by the Dyson equations 
([s]) and ^ encode the renormaUzation of the charge car- 
rier Green's functions due to the interactions with pho- 
tons and phonons, i.e. generation, recombination and 
relaxation processes. The self-energies due to interac- 
tions can be evaluated from the perturbation expansion 
of the non-equilibrium statistical average defining the in- 
teracting NEGF, using either Wick's theorem or Feyn- 
man diagrams. In the first case, the (contour ordered) 
self energy S is derived from a perturbation expansion of 
the exponential in the definition of the contour ordered 
Green's function as the nonequilibrium expectation value 
of single-carrier operators, 

G„,„' (k; i, t')^-y^ Ifc {e-i ''^^'(^)c„,k Wcl,_k(^')} 

(Al) 

a = n and k = {kx,ky,kz) (bulk) or a = i,n and 
k ~ {kx,ky) = k|| (thin film). In the following, the 
Hamiltonian terms for the perturbative interaction of 
electrons with photons and phonons shall be discussed 
and used in the derivation of the corresponding self- 
energy expressions for bulk and thin-film systems. 



1. Interaction Hamiltonian 

a. Electron-photon interaction For the electron- 
photon interaction, the perturbation Hamiltonian is 
given via the linear coupling to the vector potential op- 
erator of the electromagnetic field A, 



H.-v A • p (A2) 

mo 



with p the momentum operator and 

A(r,0 =^ [Ao(A,q)S;,,q(0+AS(A,-q)6l _q(t)" 

X e^'^^ (A3) 

where eAq is the polarization of the photon with wave 
vector q and energy hiOq added to or removed from pho- 
ton mode (A, q) by the bosonic creation and annihilation 
operators 



(A5) 



and V is the absorbing volume. 

For numerical evaluation, the Hamiltonian is reformu- 
lated in a suitable representation using the field operators 
given in (10 1 (bulk) and (461 (thin film), respectively, 



ne^{t)= J d^r¥{r,t)Hej^{r,t). 



(A6) 



In the bulk case, the above expression results iiP 



W = E E E k', q, A) 

q,A n,m k.k' 



c(i)Cmk'(0 



^A,qe 



A,-q'=^ 



, (A7) 



where the matrix element for interband transitions (n ^ 
m) is obtained from a k • p - type approximatiorP^, 



M^„(k,k',q,A) 



mo 



-^o(A,q) 



X / dVV';k(i-)(e^'"'eAq-p)7/'™k'(r) 



mo 



-Ao(A,q)(5(k' + q-k)(eAq- 



(A8) 
(A9) 



with the Bloch function momentum matrix element 



^-Wnko(f)P"mk^(f) 



(AlO) 



where f2 denotes the unit-cell volume. This gives the final 
bulk expression 



^e^W = EEE-^^^™(k'q.^) 

q.A n,m k 



X c: 



,(i)Cmk(0 



^A,qe 



A,-q'=^ 



(All) 



For devices with broken translational invariance in the 
transport dimension, the representation of the electron- 
photon Hamiltonian ( A2 ) in the real-space effective mass 



basis ( 47 1 acquires the similar from 



^e,(i) = E E E-^»V(kii'q'^) 

q.A injm kn 



bx.qC 



•i 

(A12) 



where 



■^^.^(kihq, A) 



with 



mo 



Ao(A,q)Mij(52)(eAq • p„m), 
(A13) 



M,,(g,) - / dzx*(z)e"'^-x,(^) = e'^'^'h 



(A14) 



6. Electron-phonon interaction The vibrational de- 
grees of freedom of the system are described in terms of 
the coupling of the force field of the electron-ion potential 
Vei to the quantized field U of the ionic displacemenfPSl, 



Hepir,t) = ^ W(L + /«,<) •VK.ir- (L + k)], (A15) 
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with the displacement field given by the Fourier expan- 
sion 

Uo,{Lk, t) = ^Z^„^(A, Q)e'Q-(L+-) [a^^dit) + ajv,_Q(t)] , 

{a^x,y,z), (A16) 

where the ion equilibrium position is L + k, with L the 
lattice position and k the relative position of a specific 
basis atom at this lattice site, and aA,Q,aAQ ^'"^ 
bosonic creation and annihilation operators for a (bulk) 
phonon mode with polarization A and wave vector Q in 
the first Brillouin zone. The potential felt by electrons in 
heterostructure states due to coupling to bulk phonons 
can thus be written as 



AQ 



(A17) 

where r is the electron coordinate, and ?7a,q are related 
to the Fourier coefficients of the electron-ion potentiaP^. 

The effective-mass Hamilto nian for electron-phonon 
interaction is obtained from (A17I in analogy to the 



electron-photon interaction, with the bulk result 

Q,A nM 



A,-Q 
Ua.q 



(A18) 



(A19) 



with 

7W"f(Q,A) 
and the thin-film expression 

^ep(i) = E E E ■^^(Q, A)cLk„ (oa»k„ -Q„ it) 

Q,A n k||,z 



X [oA^qe" 



a 



t 

A,-Q"= 



where 



A^f (Q,A) - ^e*0- 



(A20) 



(A21) 



For the T — X intervalley scattering considered in the 
present discussion, the coupling reads 



IC^A.QP = 



(A22) 



where a labels the phonon mode, Diy is the associated de- 
formation potential and K denotes the momentum trans- 
fer required for the scattering between two valleys. 



2. Self-energy 

At this stage, any renormalizing effect of the electronic 
system on the photons and phonons is neglected, i.e. 



the coupling to the bosons corresponds to the connec- 
tion to corresponding equilibrium reservoirs. While this 
treatment is generally a good approximation in the case 
of phonons, it is valid for the coupling to the photonic 
systems only in the case of low absorption, i.e. weak 
coupling or very short absorber length. The equilibrium 
propagators for non-interacting photons and phonons in 
isotropic media have the common form (a = 7,p) 



D'^^'^iq; E)=- 2m Nl^5{E T fu^^ 



D 



aM/A 



1 



1 



E — hujq ± if] E + huj^ ± if] 



(A23) 
(A24) 



In the above expressions, N'^^ denotes the occupation of 
the respective equilibrium boson modes, with the phonon 
occupation given by the Bose-Einstein distribution at lat- 
tice temperature T, 

^A.Q = (e^'''''-'' P = {kBT)-\ (A25) 

and the photon occupation N'l^ introduced in Sec. 
|IIB 1[ Inserting these propagators the Fock-term of 
the a general electron-boson self-energy in the first self- 
consistent Born approximation provides the steady-state 
components {a = 7,p) 



(k; E)^Y.^'^ ^) [^a" qG^ (k; T ^A,q) 

+ (Af,"q + l)G^(k„;i?±fi.;;,,q)] 
xM^^"(k,i,-q,A) (A26) 



and 



(k;i?) 



dE' S> (k;ii;')-S< (kj|;£;') 



2tt 



E' - E±iT] 



(A27) 



with the same definitions of indices and momentum as in 



(Al). Using a constant mode specific coupling strength. 



the intervalley bulk (thin-film) electron-phonon scatter- 
ing self-energy is further simplified to the diagonal form 
used in the simulations. 



<T b'^b 



2pncA 



+ iNP + l)Gf^.^,,{k;EThn,) 
where b,b' G {F^X}. 



(A28) 
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